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Abstract 

We study global features of proton-neutron (p-n) interactions and symmetry en- 
ergy over a wide range of nuclei, using a schematic model interaction with four 
forces proposed recently. Calculations are performed by the BCS approximation in 
A^, Z=20-50 and N, Z=50-82 regions. The experimental double differences of bind- 
ing energies and symmetry energy are reproduced quite well. It is shown that the 
isoscalar p-n interactions with all J are indispensable for explaining the binding 
energies of not only N ^ Z but also N > Z nuclei in the A =40-160 region. 
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In recent years, the study of proton-neutron {p-n) interactions has been a hot 
topic in the proton rich side of nuchde table. The p-n interactions are ex- 
pected to become strong in N Z nuclei where valence protons and neutrons 
in the same orbits have large overlaps of the wavef unctions, and to play very 
important roles in the structure of these nuclei [1] . The double differences of 
binding energies, which were investigated for a different purpose before [2-4], 
have been recognized as useful measures to study the p-n interactions [5-7]. 
Their observed values display characteristic behaviors near N = Z and reveal 
properties of the p-n interactions. We have analyzed what p-n interactions 
cause the characteristic behaviors in a previous paper [7]. One of our conclu- 
sions is that the double differences of binding energies trace a curve 40/^4 on 
the average and the average value is attributed to the isoscalar p-n interaction, 
almost to the J independent isoscalar p-n force. The investigation, however, 
is carried out in the N Z nuclei where a large j subshell dominates, i.e., 
f-j/2 and shell nuclei. On the other hand, the curve AQ/A observed in the 
double differences of binding energies continues till heavy nuclei. This suggests 
that considerably strong isoscalar p-n interaction persists over a wide range 
of nuclei with N > Z. 

The competition of the isoscalar p-n interaction with the isovector one is ob- 
served in the energy spectrum as well as the binding energy [8-17]. It was dis- 
cussed that the degeneracy of T = and T = 1 states in odd-odd N = Z nuclei 
are explained by a near balance of the isovector pairing and symmetry energy 
[14]. The experimental data indicate that the symmetry energy accompanied 
by the so-called Wigner energy is proportional to T(T-\-l) with a common co- 
efficient a(A),i.e., Esym + Ew = a{A)T (T+l)/ A where T = \Tz\ = \N-Z\/2. 
It has been shown that the Wigner energy is originated in the isoscalar p-n in- 
teraction [18,19] but cannot be solely explained in terms of the J = 1 isoscalar 
p-n interaction [15,19]. These investigations have suggested the importance of 
isoscalar p-n interactions with various J in the symmetry energy and Wigner 
energy. This is confirmed in the single j shell model calculation for /7/2 shell 
nuclei in our previous paper [7]. Our model reproduces very well the empiri- 
cal symmetry energy coefficient a{A). The result indicates that the symmetry 
energy coefficient is explained by the J independent isoscalar p-n force and 
J — isovector pairing force. In this letter, we extend our discussion to a 
wide range of nuclei. The empirical formula of Seeger [20,21] or Dufio and 
Zucker [22] for a{A) has a mass dependence (1 — r]A~^^'^). This A dependence 
combined with the curve 40/ A of the average double differences of binding 
energies will give a good guide to determine the A dependence of our isoscalar 
p-n force strength. 

Our previous study [7] has been carried out using a schematic interaction 
composed of four forces which is extended from the pairing plus QQ force. 
This schematic interaction reproduces well the low-lying states and binding 
energies oi N ^ Z nuclei where a large j subshell dominates, and also var- 
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ious nuclear properties, such as the double differences of binding energies, 
symmetry energy, Wigner energy, odd-even mass difference and two-proton 
separation energy. The p-n interactions play leading roles in these properties, 
especially the J independent isoscalar p-n force is essential for the properties 
which we consider here. The schematic interaction with four force parameters 
is easy to extend the model space. In this letter, we investigate systematics 
of the double differences of binding energies and symmetry energy in both of 
N ^ Z and N > Z nuclei with N, Z =20-50 and N, Z =50-80. We perform 
this many j shell calculation in the BCS approximation. 

The Hamiltonian we proposed in Ref. [17] has an isospin-invariant form as 
follows: 



H = Hsp + Hp^ + Hqq + Hp^ + °, (1) 
^sp = - A)4^c«p, (2) 

ap 

Hpj = -^gj Y.Y.PA(ib)pj{cd)A\Mi^{ab)AjMin{cd), (3) 

MKa<bc<d 

Hqq = —X mm q{ah)q{cd) : BlMp{ah)B2Mp'{cd) : (4) 

^ M abp cdp' 

= E E AlMooi<^b)AjMoo{ab), (5) 

a<bJM 



with 



^jMr«(a&) = B\j^p{ab) = [cipCbp]jM- (6) 

Here po{ab) = ^ {2ja + l)Sab and p2{ab)=q{ab) /b'^={a\\r^Y2/b'^\\b) / y/E where 6^ 
is the harmonic-oscillator range parameter. This Hamiltonian has only four 
parameters, go, g2, X, and k^. 

We consider the model space (I/7/2! 2p3/2, if 5/2, 2pi/2, 1(79/2) for N, Z = 20—50 
region and (2^5/2, l5'7/2i 3si/2! 2(i3/2, l/iii/2) foi' N, Z = 50—82 region. The low- 
est single-particle energy in each region, is determined as ej^^^ — ^ = -B(^^Ca) — 
5('i0Ca) = -8.3633 MeV and e^^^^ - A = BC^'^Sn) - B(^^'^Sn) = -10.265 MeV 
from the experimental binding energies. The other single particle energies are 
chosen so as to be the same level spacings as those used by Kisslinger and 
Sorensen [23]. The pairing force strength Qq and the quarupole force strength 
X are ajusted to fit approximately the experimental odd-even mass differences 
and the lowest 7 = 2"*" energies in quasiparticle Tamm-Dancoff approximation, 
respectively. 
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go^2AA-\ x = 350A-5/V^ for AT, Z = 20 - 50 region, (7) 
go^2AA-\ x = 450A-^/V^ for AT, Z = 50 - 82 region. (8) 



The quadrupole pairing force strength g2 can be determined by the same 
relation g2 — 0.2go as that used by Hara et al. [24]. We shall fix the strength 
of the isoscalar p-n force (5) in the following discussion. 

Let us first analyze the experimental values of the mth double difference of 
binding energies defined by the following equation [7] , 

5\/(™)(Z,iV) = (5(™)5(Z,iV), (9) 

where B[Z, N) is the nuclear binding energy and the operator 5^™^ is defined 
as 



N) = -^[f{Z, N) - f{Z, N-m) 

-f{Z -m,N) + f{Z -m,N -m)]. (10) 

Figure 1(a) shows the plot of 5V^^\Z,N) as a function of ^ = + Z for 
nuclei in the mass region ^4=16-210. We can sec two separate groups in the 
figure, namely, one is for the even A nuclei (dots) and the other is for the odd 
A nuclei (crosses). The patterns of dots and crosses show shell effect at Z or 
N =28, 40, 50, 82 and are symmetric with respect to the average curve 40/A. 
The double difference of binding energies 5V^^\Z,N) can be approximately 
written as 



5y«(z,Ar) = ^ + (-i)V, (11) 

where /' denotes a deviation from the average 4Q/A. The large values of 
5V^^\Z,N) for even-A nuclei (circled dots) ai N = Z are striking below 
A = 80. The data of another double difference of binding energies dV^'^^Z, N) 
are plotted in Fig. 1(b) as a function oi A — N + Z. Our definition of SV^"^^ 
has the opposite sign to Brenner's one [6]. The SV^'^'' values show a different 
behavior from the 6V^^^ ones. We see large scatters of dots and crosses for 
A < 80, which show of iV ^ Z nuclei. The values of N = Z are 

especially large as that of 6V^^\ If we neglect the scattered dots and crosses 
m N !v Z nuclei, SV^'^^ varies rather smoothly. This smooth trend is clear for 
A > 80 and continues up to heavy nuclei. This is due to the fact that there 
are no stable N ^ Z nuclei with A > 80. The systematic decrease of 6V'^'^^ 
with increasing mass A can be traced by the same average curve AO /A as that 
of dV^^\ The deviations from the curve 40/^4 are small and shell effect disap- 
pears for 6V^'^\ In Ref. [7], we analyzed 6V^^^ and 5V^'^^ for probing the p-n 
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interactions in N ^ Z nuclei of 5^9/2 and /7/2 subshells. The analysis reveals 
different contributions of isoscalar and isovector p-n interactions to the double 
differences of binding energies. The average curve AO /A of dV^^^ and dV^"^^ is 
attributed almost to the isoscalar p-n force introduced, while the p-n 

part of Pq + QQ + P2 force produces the deviation from the curve AO /A in 

N). 

bmce comes from the J independent isoscalar p-n force H^^, the average 
curve AO/ A of 5V^'^^ reflect the A dependence of the force strength k^. The 
interaction energy of H^"^^ in the state with the total valence nucleon number 
n and isospin T is given by 

1 71 71 

K:' - {H:f) = + 1) - r(r+ 1)}. (12) 

Table 1 The coefficient of empirical symmetry energy. 



Osym (MeV) 


References 


28.1 
23.3 
25.0 

20.65 (1.0 - 2.32/y4V3) 

30.59 (1.0 - 1.76/^V3) 

33.60 (1.0 - 1.52/^V3) 


Ring and Schuck [25] 
Droplet model by Myers [26] 
Bohr and Mottelson [27] 
Green [28], Seeger [20] 
Seeger [21] 
J.Duflo and A.RZucker [22] 



If we calculate the double difference of the interaction energies f/^^ ^ (which 
we denote by S^'^^U^^^) by substituting t/^j7° for / in Eq. (10), we can expect 
5(2)f;r=o ^ ^y(2) discusscd in Ref. [7] that the average curve 40 /A of 
5y(2^ suggests the force strength — 80/A in the lowest order. The simple 
form = 80/^4, however, cannot well reproduce the curve 40/^4 for S^'^^U^^^ 
as shown in Fig. 2. The disagreement demands a higher order correction. 
According to our analysis [7], the symmetry energy is also governed by the 
isoscalar p-n force Let us look at the symmetry energy coefficient agym ~ 

0(^4)74 in order to get a hint on the A dependence of Table 1 hsts the 
values of ttsym proposed by various authors. The elaborate coefficient in Rcfs. 
[20-22,28] has a mass dependence as the form (1 — r]/A^^'^) (the lowest A 
dependence l/A of the symmetry energy is not included in the definition of 
<^sym)- The following form of /c'^, therefore, deserves to be tested: 



k'^^{l.O-r^/A'/% 



(13) 



where ( and r] are parameters. Figure 2 shows the values of 5*^^^f/^j7° calculated 
with the use of the force strength (13). The calculated values finely 

reproduce the curve 40/^4 when we choose the parameters ( = 224.0, r] = 2.2 
for the TV, Z = 20 - 50 region and C = 224.0, 77 = 3.0 for the A^, Z = 50 - 82 
region. 

The four parameters {go,Xj 92, k^) of our Hamiltonian were fixed by Eqs. (7), 
(8) and (13). We calculated the binding energies E of A = 40 — 140 nuclei 
within the BCS approximation, i.e., the expectation values of the Hamiltonian 
with respect to the BCS states for the A systems. (We did not calculate the 
binding energies near N = Z nuclei and near the beginning and end of each 
shell region in order to avoid the error of the BCS approximation.) The double 
differences of binding energies calculated from the approximate ground-state 
energies E are shown in Figs. 3(a) and 3(b). The calculated values are nicely 
around the curve AO /A. Especially, 6'^'^^E reproduces quite well the observed 
variation of S^'^^V in Fig. 3(b). We can say that our schematic Hamiltonian can 
qualitatively explain the behavior of the double differences of binding energies 
also in many-j shell case. 

The double differences of binding energies, 5^^^V and S^^^V, are considered 
as good measures of the effective p-n interactions. In the previous paper [7], 
we investigated from what p-n inteactions the characteristic behavior of d^^^V 
and come, in terms of the single-j shell model for N ^ Z nuclei. Let 

us develope the analysis in the present many-j' shell calculation. We consider 
the contributions of respective forces to the ground-state energy in even-even, 
even-odd, odd-even and odd-odd nuclei, respectively, 

Eie,e) = {0\H\0) = Uo + U:=', (14) 

Eie, o) = {n\H\n) = f/o + A„ + [/;=°, (15) 

E{o, e) = {p\H\p) = Uo + Ap + U^f, (16) 

E{o, o) = {pn\H\pn) = [/q + + A„ + e^n + Kf. (17) 

In the BCS approximation, the ground states of the four types of nuclei given 
above are defined by the quasiparticle vacuum, one quasiparticle and p-n two 
quasiparticle states, 

|0), |n) = 4|0), |p) = 4|0), M = aX|0), (18) 

We denote the ground-state energy of the even-even system |0) by Uq, the 
neutron (proton) gap by A„ (Ap) and the interaction energy of the residual 
force {QQ + P2) by 

Uo = {0\H,, + Hp, + Hqq + Hp,\0), (19) 
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/:^p = {p\Hsp + Hp, + HQQ + Hp,\p)-U^ = |n) or (20) 

epn = {pn\HQQ + Hp^ \pn) . (21) 

It is easy to calculate d'^^^E and from Eqs. (14-17). We have seen in 

Fig. 2 that the value of S^'^^U^^^ yields the curve 40/A, when we use the 
parameters {(, rj) fixed above. We therefore obtain the relation S^^^E(e,e)— 
S^^^E{o, o)=AO/A - €pn for even A nuclei, ^(^^^(e, o)=S^^^E{o, e) =AO/A + €pn 
for odd A nuclei and S^'^^E = AO /A both for even and odd nuclei. The char- 
acteristic behaviors of d^^^V and ^'-^-'V" are clearly understood. Their average 
values AO /A are attributed to the J independent isoscalar p-n force also in 
the present many-j shell calculation. The deviations from the average curve 
AO /A in are due to the residual p-n interaction between the last proton 

and neutron in odd-odd nuclei caused by the QQ + P2 force (the QQ force 
mainly contributes to it). It is notable that the residual p-n interaction e^^ 
is attractive. This and the relation E{o, o) — E{e, e)=Ap + A„ + ep„ obtained 
from Eqs. (14) and (17) explain the fact that the binding energy difference 
B{o, o) — B{e, e) is systematically smaller than the sum of gaps Ap + A„ (Bohr 
and Mottelson [27] suggested an attractive residual p-n interaction about 20/^4 
between the last odd proton and neutron in odd-odd nuclei). 

Let us now discuss the symmetry energy Egy^ and the Wigner energy E\y. 
Wc calculated Egy^ and Eyy from the ground-state energies in the BCS ap- 
proximation. Figure 4 shows the calculated and experimental symmetry cnrgy 
coefficient a{A) in the expression Esym + Ew = a{A)T{T + I) /A. (The A de- 
pendence of a{A) = Attsym extracted by Seeger [21] and Duflo et al. [22] are 
also shown in Fig. 4.) The figure is plotted separately for "diagonal" and "off- 
diagonal" regions where neutrons and protons are in the same major shell for 
the former and in the different major shells for the latter. We estimated the 
Coulomb-energy-corrected binding energies B* = B{exp) + i?coui(cal) follow- 
ing Comay and Janecke [31] and extracted the symmetry energy coefficients 
by the treatment of Janecke and Comay [32] . The agreement between theory 
and experiment is very well for the diagonal region of A^, Z = 20 — 50, and 
is qualitatively good for the off-diagonal region of A^, Z = 20 — 50 and the 
diagonal region of N, Z = 50 — 80 where the calculated values are a little 
smaller than the experimental ones. The coefficient a{A) in the off-diagonal 
region is approximately constant for the heavier nuclei as expected from the 
simple liquid-drop model. They are larger than those of the diagonal region, as 
pointed out by Janecke and Comay [32] . In our calculation, the single-particle 
energy Esp = {0\Hsp\0) influences the symmetry energy coefficient because of 
a large energy gap between the major shells. We can also see the shell-closure 
effect on the symmetry energy as well as the double difference of binding 
energies SV^^\ 

We analyze the symmetry energy in the same way as the single-j shell model 
calculation [7]. If the single-particle energies are roughly degenerate, i.e., 



7 



e, and Hqq + i7p„ can be neglected, the Hamiltonian Hsp + Hp^ + H^^^ has 
the SO (5) symmetry. Then, the groimd-state energy is approximated in terms 
of the total valence nucleon number n and the total isospin T as follows: 

E^{H,, + Hp,) + U;f 

= (e - X)n - -{gonin - + + 1)} 

+\{go + k')T{T+l). (22) 

The term proportional to T{T + 1) in the ground-state energy (22) just cor- 
responds to Egym + Ew- The approximate symmetry energy coefficient a{A) 
is given by 

~a{A) = ]^{go + k')A. (23) 

When we use the parameters go = 2A/A and A:° = 224(1.0 -2. 2/A^/''^)/A in the 
fpg region, the symmetry energy coefficient in the above SO (5) limit becomes 
a^yin = a{A)/A = 31(1.0 - 1.987/^1/3) from Eq. (23). This agrees considerably 
well with the empirical symmetry energy 30.59(1.0 - 1.76/AV3) extracted by 
Seeger [21]. We can, therefore, say that the symmetry energy and Wigner 
energy are originated nearly in the J independent isoscalar p-n force and the 
J = isovector pairing force. The two contributions to Ssym are, for instance, 
12.3 MeV and 3.0 MeV respectively in the A = 60 system. The isoscalar p-n 
force H^^^ plays a significant role in reproducing the experimental symmetry 
energy and the Wigner energy. 

In conclusion, the systematic behaviors of the double differences of binding 
energies, symmetry energy and Wigner energy are reproduced well by the 
schematic interaction Pq + QQ -|- P2 + H^7^ within the BCS approximation, 
over the wide range of nuclei. It was shown that the J independent isoscalar 
p-n force is crucial for nuclear properties not only in N Z nuclei but also in 
N > Z nuclei far from the N = Z line. The analysis also tells that the residual 
interaction between the last proton and neutron in odd-odd nuclei consists in 
that of the QQ force. We estimated the effect of the ground-state correlation 
in the framework of of the Random-Phase approximation. The calculations, 
however, did not change the conclusions mentioned above. 
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Figure captions 



Fig. 1 The double differences of binding energies derived from the experi- 
mental masses in the region A — 16 — 210: (a) 6V'^^^{Z, N) as a function of 
A^N + Z; (h) SV^'^^Z, AT) as a function ofA = N + Z. The dots stand 
for even A nuclei, and the crosses for odd A nuclei. The curve 40 /A is also 
plotted. 

Fig. 2 The calculated values of 6^^'>U;=^{Z, N) as a function oi A = N + Z 
in A^, Z=20-50 and A^, Z=50-82 regions. 

Fig. 3 The calculated values of S^^^E{Z, N) and S'^'^^E{Z, N) as a function of 
A = N + Z in N, Z=20-50 and N, Z=50-82 regions. 

Fig. 4 The symmetry energy coefficient in A^, Z=20-50 and A^, Z=50-82 re- 
gions. The calculated and experimental values are denoted by the open 
and filled circles, respectively. The full hues shows the empirical formula by 
Seeger [21], and the dotted line the empirical formula by Duflo and Zucker 
[22]. 
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